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ABSTRACT 

Glitches detected in pulsar timing observations at the Nanshan radio telescope of 
Urumqi Observatory between 2002 July and 2008 December are presented. In total, 
29 glitches were detected in 19 young pulsars, with this being the first detection of 
a glitch in 12 of these pulsars. Fractional glitch amplitudes range from a few parts 
in 10"" to 3.9 X 10"^. Three "slow" glitches are identified in PSRs J0631+1036, 
B1822— 09 and B1907+10. Post-glitch recoveries differ greatly from pulsar to pulsar 
and for different glitches in the same pulsar. Three small glitches in PSRs B0402+61, 
B0525+21 and J1853-I-0545 show evidence for normal post-glitch recovery, but for 
PSRs B0144-f 59 and B2224-I-65 the spin frequency v continually increases relative 
to the pre-glitch solution for hundreds of days after their small glitches. Most large 
glitches show some evidence for exponential post-glitch recovery on timescales of 100 
- 1000 days, but in some cases, e.g., PSR B1758— 23, there is little or no recovery, that 
is, no detectable increase in |z>| at the time of the glitch. Beside exponential recoveries, 
permanent increases in slowdown rate are seen for the two large glitches in PSRs 
B1800— 21 and B1823— 13. These and several other pulsars also show a linear increase 
in i/ following the partial exponential recovery, which is similar to the Vela pulsar 
post-glitch behaviour. However, the rate of increase, i>, is an order of magnitude less 
than for the Vela pulsar. We present improved positions for 14 of the glitching pulsars. 
Analysis of the whole sample of known glitches show that fractional glitch amplitudes 
are correlated with characteristic age with a peak at about 10^ years, but there is a 
spread of two or three orders of magnitude at all ages. Glitch activity is positively 
correlated with spin-down rate, again with a wide spread of values. For some (but not 
all) pulsars there is a correlation between glitch amplitude and the duration of the 
following inter-glitch interval. In no case is there a correlation of glitch amplitude with 
the duration of the preceding inter-glitch interval. 
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1 INTRODUCTION 

After allowing for the accurate determination of the pul- 
sar spin-down, pulse arrival-time measurements display two 
types of irregularities in the pulsar rotation rate: timing 
noise and glitches. Timing noise is a continuous wandering of 
the pulsar rotation rate, on time-scales of days, months and 
years. In most young pulsars, timing nois e is generally very 
"red" , i.e., stronger at lower frequencies (|Cordes &: Downa 
[l98^, and the second time-derivative of the pulse frequency 
is commonly used as an indicator of the level of noise. 

Glitches are characterized by a sudden increase in the 
pulsar rotation rate, often followed by a (partial) recov- 
ery toward the pre-glitch status. Glitches were first ob- 
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serve d in the Vela pulsar ijRadhakrishnan fc Manchester! 
1 19691 : iReichlev fc Downs! ll969D and it was soon realized 
that they could be an important diagnos tic for studyin g 
neutron-star interiors (Rudcrman 1969; Ba vm et al.!ll969l ). 
It is widely believed that these events are caused by either 
sudden and irregular transfer of angular momentum from a 
faster rotating interior super-fluid to the solid crust of the 
neutron star, or quakes in a solid component of the star, 
or even both. The result is a sudden increase in the rota- 
tional frequency u of the pulsar with a relative magnitude 
normally in the range 10"^" < Aug/u < 5 X 10"*'. Large 
glitches {Aug/u > 10^^) are usually accompanied by jumps 
in the spin-down rate {Avgjv ~ 10~* to 10~^). A char- 
acteristic feature of many glitches is a relaxation after the 
frequency jump, which may occur over periods of minutes to 
years. In some cases, this relaxation is approximately expo- 
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nential toward the extrapolated pre-glitch state, but many 
different forms of relaxation are o bserved, some with mul- 
tiple timescales. iLvne et alJ (2000) presented evidence that 
the degree of relaxation in different pulsars is related to their 
spin-down rate 

Although glitches are a rather rare phenom- 
in creasing number of these events have been 



ena, an 



reported JWang et all |2000l : iKrawczvk l2003l : 



Ijanssen fc Stapper^ 2006f K p ermitting 



more compre- 



hensive statistical studies fe.g.. [Lvne et al.ll2000l ). To date 
about 170 glitches have been observed, in 51 of the over 
1800 known pulsar^. The most frequent glitches have 
been detected in PSRs B053H~21 (Crab), J0537-6910, 
B0833-45 (Vela), J1341-6220 and B1737-30, which have 
characteristic ages in the range ~ 10^ to 10^ yr. PSR 
J0537— 6910 has glitched 23 times in seven years, with 
most of the glitches being relatively small in ^Vgjv terms 
(typically a few parts in 10~^). Their amplitude shows a 
strong correlation with t he interval to the following glitch 
l|Middleditch et al.ll2006l '). Glitches in the Crab pulsar are 
even smaller, with A.Ug/v typically ~ 10~^ (jLvne et alJ 
I1993D . There is normally a rapid exponential decay over 
several days and a persistent increase in slow-down rate 
l|Lvne et alJ Il993l ). The Vela pulsar is the prototypical 
pulsar showing large glitches, with /S.Vg/v generally > 10~® 
although some small and intermediate-sized glitches have 
been observed. The post-glitch decay is initially exponen- 
tial, with evidence for several timescales ranging from 1 
minute to several hundred days, but subsequently linea r 
till the next glitch (|Dodson et al.ll2007l : [Lvne et ai.lll996l ). 
PSR B1737— 30 is a frequently glitching pulsar which shows 
a large range of glitch sizes and a high value of the glitch 
activity parameter 
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where T is the total data span (iMcKenna fc Lvn3 Il990l : 
IShemar fc Lto3 1 19961 : IZou et al.l |2008| ). Significant post- 
glitch decays are observed for the larger glitches. 

Clearly, the glitch phenomenon in pulsars is complex. 
Further observations of a larger sample of pulsars will help to 
quantify and understand processes responsible for the very 
different glitch and post-glitch behaviours observed both be- 
tween different pulsars and in different glitches in the same 
pulsar. 

Nearly 280 pulsars are being monitored regularly at 
Urumqi Observatory in a timing program which commenced 
in 2000. The aims of the project are to make precise mea- 
surements of the position and timing behaviour of pulsars, 
to provide ephem eris for observatio ns at radio and other 
wavelengths (e.g.. [Smith et al]l2008l ) and to investigate the 
interiors of neutron stars by monitoring irregularities in their 
rotation behaviour. In this paper we describe the glitches de- 
tected in our observations, together with updated descrip- 
tions of the relaxations for some previously reported glitches. 



2 OBSERVATIONS AND ANALYSIS 

Timing observations of pulsars using the Urumqi Obser- 
vatory 25-m radio telescope at Nanshan started in 2000 
January with a dual-channel room-temperature receiver. A 
dual-channel cryogenic receiver system has been used at 
the central observing frequency of 1540 MHz since 2002 
July. Approximately 280 pulsars are observed regularly three 
times each month and results to 2008 December are reported 
here. The two polarizations, each of bandwidth 320 MHz, are 
sent to a filter-bank consisting of 2 x 128 channels of width 
2.5 MHz. The data are digitized to one-bit pr ecision, and 
the i ntegration time is typically 4 to 16 min (|Wang et al.l 
I2OOII ). 

Off-line data reduction was performed in the follow- 
ing steps. The data for each pulsar were dedispersed and 
summed to produce a total intensity profile centered at 
1540 MHz. Local arrival times were determined by correlat- 
ing the data with standard pulse profiles; the pulse times of 
arrival (TOAs) normally correspond to the peak of the main 
pulse. TOAs at the Solar-system barycenter were obtained 
using the standard timing program TEMPCEI with the Jet 
Propulsion Laboratories DE405 ephemeris l|Standishlll99g '). 
A model for the pulsar timing was fitted to these barycentric 
TOAs, weighted by the inverse square of their uncertainty. 
Errors in the fitted parameters are taken to be twice the 
standard errors obtained from TEMPO. The basic timing 
model for the barycentric pulse phase, <j), as a function of 
time t, is 

4>(t) ^<t>o^ v{t - to) + ^i>{t - tof + ii>(f - ^o)^ (2) 

where (fiQ is the phase at time to and v, i>, i) represent 
the pulse frequency, frequency derivative and frequency sec- 
ond derivative respectively. Pulsar positions were derived 
from long inter-glitch intervals with the residuals being 
"whitened" by fitting of higher-order frequency derivatives 
if necessary. 

Glitches can be described as combinations of steps in 
v and v, of which parts may decay exponentially on vari- 
ous timescales. We assume a glitch model for the observed 
phenomenon as below: 

v{t) = vo{t) + Ai/p + A!>pt + Ai^de"*/^'' 



(3) 
(4) 



where Avp and Avp are, respectively, the permanent changes 
in frequency and frequency derivative relative to the pre- 
glitch values and Aud is the amplitude of the exponential 
recovery with a decay time constant of Td- We note that for 
glitches the jump in frequency 1/ is positive and the jump in 
frequency derivative z> is usually negative. Since u is nega- 
tive, a negative jump corresponds to an increase in magni- 
tude of the spin-down rate. The total frequency change at 
the glitch is 

Aug^Aup + Aud. (5) 
The degree of recovery is often described by the parameter 



1 See the ATNF Pulsar Catalogue l|Manchester et al.l l2005h 
Glitch table at www.atnf.csiro.au/research/pulsar/psrcat 



Q = AVdjAVg 



See |http:// www, atnf . csiro.au/ research / pulsar / tempo | 
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Table 1. Parameters for the pulsars with observed glitches 



Pulsar Name 


RA 


Dec 


Epoch 


P 


P 


Age 




References 






(h m s) 


(d m s ) 


MJD 


(s) 


(10-15) 


(kyr) 


(IQi^G) 


for position 


J0147+5922 


B0144+59 


01:47:44.6457(17) 


+59:22:03.293(15) 


53160 


0.1963 


0.256 


12100.0 


0.23 


this work 


J0406+6138 


B0402+61 


04:06:30.082(3) 


+61:38:41.04(3) 


53784 


0.5945 


5.576 


1690.0 


1.84 


this work 


J0528+2200 


B0525+21 


05:28:52.264(9) 


+22:00:04(2) 


53100 


3.7455 


40.036 


1480.0 


12.40 


this work 


J0631+1036 




06:31:27.524(4) 


+ 10:37:02.5(3) 


53850 


0.2877 


104.683 


43.6 


5.55 


this work 


J1705-3423 




1 7 05-42 363('3'1 


— 34-2'^-45 17C121 


53830 


0.2554 


1.075 


3760.0 


0.53 


this work 


J1740-3015 


B1737-30 


17:40:33.82(1) 


-30:15:43.5(2) 


52200 


0.6068 


466.358 


20.6 


17.00 


1 


J1751-3323 




17:51:32.725(11) 


-33:23:39.6(9) 


53750 


0.5482 


8.897 


976.0 


2.23 


this work 


J1801-2304 


B1758-23 


18:01:19.829(3) 


-23:04:44.2(2) 


50809 


0.4158 


112.880 


58.4 


6.93 


2 


J1803-2137 


B1800-21 


18:03:51.411(1) 


-21:37:07.35(1) 


51544 


0.1336 


134.100 


15.8 


4.28 


3 


J1812-1718 


B1809-173 


18:12:07.208(8) 


-17:18:29.5(11) 


49612 


1.2053 


19.770 


1000.0 


4.85 


4 


J1818-1422 


B1815-14 


18:18:23.7660(17) 


-14:22:36.71(16) 


54000 


0.2914 


2.038 


2270.0 


0.78 


this work 


J1825-0935 


B1822-09 


18:25:30.629(6) 


-09:35:22.3(3) 


53300 


0.7689 


52.285 


233.0 


6.42 


this work 


J1826-1334 


B1823-13 


18:26:13.175(3) 


-13:34:46.8(1) 


52400 


0.1015 


75.061 


21.4 


2.79 


this work 


J1841-0425 


B1838-04 


18:41:05:663(1) 


-04:25:19.50(7) 


53900 


0.1861 


6.391 


461.0 


1.10 


this work 


J1853+0545 




18:53:58.411(1) 


+05:45:55.26(3) 


53882 


0.1264 


0.611 


3280.0 


0.28 


this work 


J1902+0615 


B1900+06 


19:02:50.277(2) 


+06:16:33.41(6) 


53390 


0.6735 


7.706 


1380.0 


2.31 


this work 


J1909+1102 


B1907+10 


19:09:48.694(2) 


+ 11:02:03.35(3) 


49912 


0.2836 


2.639 


1700.0 


0.88 


4 


J1915+1009 


B1913+10 


19:15:29.984(2) 


+10:09:43.67(4) 


53300 


0.4045 


15.251 


420.0 


2.51 


this work 


J2225+6535 


B2224+65 


22:25:52.721(7) 


+65:35:35.58(4) 


53880 


0.6825 


9.659 


1130.0 


2.58 


this work 



llFomalont et al.l lll997h: 2. l^aUetal] 1IT99I); 3.lBrisken et al.l (120061): 4.lHobbs etHI ll2004h. 



Because of the decaying component, the instantaneous 
change in v at the ghtch differs from Ai>p 

Ai>g = Ai>p + A!>d = A!>p - Qd^vg/n (7) 

Solutions before and after the ghtch are used to estimate 
the steps in frequency and its derivative. If the gap be- 
tween observations around the glitch is not too large, the 
glitch epoch can be estimated more accurately by requiring 
a phase-connected solution over the gap in the TEMPO fit. 
However if the gap is too large, rotations can be missed and 
the glitch epoch is estimated as the mid-point of the gap 
between observations around the glitch. 



3 RESULTS 

We present the results of glitch analysis using Nanshan tim- 
ing data from July 2002 to Dec 2008. Table [T] lists the pa- 
rameters for the 19 pulsars for which glitches were observed. 
The first two columns give the pulsar names based on J2000 
and B1950 coordinates and the third and fourth columns are 
their J2000 positions. The next four columns give the pulsar 
period and period derivative, characteristic age Tc — P/{2P) 

and surface dipole magnetic field Bs = 3.2 x lO^' \/ PP G. 
Most of the positions quoted in Table [T] result from our work 
except for five pulsars where the referenced positions are 
more accurate. Estimated uncertainties in the last quoted 
digit, given in brackets, are 2a from the TEMPO fit in this 
and following tables. References for the position are given 
in the final column. The data spans used to determine the 
positions are the largest available which are not too strongly 
affected by post-glitch recovery. Typical spans are 3 - 5 yr, 
with a minimum of 2.3 yr and a maximum of 7.5 yr. 



As mentioned in Section[2]the data sets used for position 
fitting are whitened by fitting for higher order derivatives of 
frequency. In most cases, five or less derivatives are sufficient 
to whiten the data, giving rms residuals of 100 - 700 /^s. For 
PSR J1751-3323, nine derivat ives were r e quired and for PSR 
B1815-14, seven derivatives. IZou et al] l|2005h investigated 
the effect of timing noise on the determination of pulsar 
positions, showing that an uncertainty of twice the formal 
TEMPO error is generally realistic for data spans of the 
order of 5 yr. 

In total, 29 glitches were detected in these pulsars. The 
rotational parameters from independent fits to the pre-glitch 
and post-glitch data are given in Table [21 Except for short 
sections of data, the fits include a u term, which is gener- 
ally dominated by timing noise or recovery from a previous 
glitch. 

Glitch parameters are given in TableO in which the sec- 
ond column gives the glitch epoch, and the next two columns 
give the fractional changes in v and u at the glitch deter- 
mined by extrapolating the pre- and post-glitch solutions to 
the glitch epoch. Quoted uncertainties include a contribu- 
tion from the uncertainty in glitch epoch. Results of directly 
fitting for the glitch parameters in TEMPO are given in the 
remaining columns. Decay terms were only fitted when sig- 
nificant decay was observed. 

The rotational behaviour of these pulsars are shown in 
Fig. [Tlto rill Values of u and were obtained from indepen- 
dent fits of Equation [2] (omitting the i> term) to data spans 
ranging from 50 d to 200 d. The residual frequencies have 
been obtained at the various epochs by the subtraction of 
the pre-glitch models. We discuss the results for each pulsar 
in the following Sections. 
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Table 2. Pre- and post-glitch timing solutions 



Pulsar Name 




i) 


V 


Epoch 


MJD Range 


No. of 


RMS 




(s-i) 


(10-12 s-2) 


(10-2* s-^) 


(MJD) 






TOAs 


(/is) 


"r~> f\t A A 1 ^ f\ 

B0144-|-o9 


5.093689334287(4) 


—0.0066633(2) 




53160.0 


CO/1 QC 


-53662 


90 


101 




5.093688844679(4) 


-0.0066619(2) 




54011.0 


53686- 


-54831 


107 


145 




1.08187172559(1 ) 


—0.015757(2) 




52745.0 






57 


395 




1.681870311647(2) 


-0.01576474(8) 




53784.0 


53050- 


■54830 


189 


541 




0.26d984o2o978(2) 


—0.0028537(3) 




51864.0 


51546- 


-52271 


85 


575 




0.266984598507(1) 


—0.00285481(3) 




52800.0 


52287— 


-53975 


249 


834 




0.266984253314(2) 


-0.0028550(1) 




54200.0 


53982- 


-54831 


182 


817 


J0631-t-1036 


3.47486317481(6) 


—1.263638(6) 


30(2) 


52657.0 


52473- 


-52842 


40 


274 




3.47482575938(6) 


—1.265058(9) 


106(1) 


53000.0 


ozttoo— 


-Ooziy 


37 


685 




3.474732970746(8) 


-1.263032(1) 


5.47(8) 


53850.0 


53238- 


-54100 


175 


1454 




3.47467842199(3) 


-1.263630(2) 


15.7(6) 


54350.0 


54129- 


-54630 


106 


5003 




3.4746348982(1) 


-1.26374(2) 


166(12) 


54750.0 


54634- 


-54831 


57 


315 


J1705-3423 


3.915021810541(8) 


-0.0164948(2) 


— 


53800.0 


52486- 


-54384 


222 


670 




3.9150208150(1) 


-0.016495(9) 


- 


54500.0 


54394- 


-54832 


44 


590 


B1737-30 


1.647839443329(3) 


-1.2659499(3) 


13.16(9) 


54000.0 


53101- 


-54394 


115 


1962 




1.64778155382(3) 


-1.26561(1) 


59(3) 


54530.0 


54463- 


-54685 


19 


314 




1.6477542263(1) 


-1.26557(3) 


- 


54780.0 


54695- 


-54808 


7 


338 


J1751-3323 


1.82406604540(3) 


-0.029722(4) 


_ 


52758.0 


52496- 


-52994 


35 


1121 




1.82406350461(2) 


-0.0297075(4) 


-0.58(4) 


53750.0 


53007- 


-54384 


75 


3667 




1.8240615852(2) 


-0.02938(2) 




54500.0 


54463- 


-54714 


18 


1118 


B1758-23 


2.40490661705(4) 


-0.65316(2) 


3(3) 


52900.0 


52494- 


-53292 


67 


3079 




2.40484572830(1) 


-0.6531075(6) 


1.12(4) 


54000.0 


53327- 


-54831 


116 


2604 


B1800-21 


7.48196609098(1) 


-7.4862079(4) 


238.57(3) 


52800.0 


51916- 


-53429 


165 


7091 




7.481282565107(9) 


-7.5200390(9) 


207.93(4) 


53900.0 


53458- 


-54832 


108 


12533 


B1809-173 


0.82961975150(3) 


-0.013100(4) 




52800.0 


52529- 


-53105 


32 


1360 




0.82961862949(2) 


-0.0131283(4) 


0.25(5) 


53800.0 


53119- 


-54351 


79 


1013 




0.82961772349(6) 


-0.01313(1) 


- 


54600.0 


54382- 


-54764 


17 


1214 


B1815-14 


3.43066049095(2) 


-0.024011(3) 


— 


51822.0 


51512- 


-52048 


69 


558 




3.430655769877(4) 


-0.0239594(2) 


-0.446(3) 


54100.0 


52081- 


-54831 


226 


2400 


B1822-09 


1.300387566396(6) 


-0.0886064(2) 


4.34(4) 


53280.0 


52811- 


-53716 


83 


1362 




1.30038260224(2) 


-0.088638(3) 


-22.0(6) 


53930.0 


53723- 


-54094 


45 


1203 




1.30038021608(2) 


-0.088780(2) 


-2.8(1) 


54262.0 


54129- 


-54831 


54 


1543 


B1823-13 


9.85460718889(1) 


-7.2729146(3) 


105.74(2) 


52400.0 


51574- 


-53199 


132 


3407 




9.85394781571(6) 


-7.265061(9) 


66(1) 


53450.0 


53254- 


-53730 


32 


389 




9.85349875132(2) 


-7.3062994(9) 


280.88(9) 


54200.0 


53739- 


-54518 


90 


14061 


B1838-04 


5.37205530549(3) 


-0.184289(2) 


1.5(2) 


53044.0 


52535- 


-53387 


57 


364 




5.37204477631(1) 


-0.1844913(5) 


-1.56(4) 


53900.0 


53429- 


-54662 


99 


556 


J1853-I-0545 


7.91137412600(2) 


-0.038303(1) 


- 


53100.0 


52493- 


-53434 


63 


446 




7.91137154739(1) 


—0.0383351(1) 


—0.19(5) 


53882.0 


53458- 


-54830 


113 


549 


B1900H-06 


1.484775165625(2) 


-0.01700231(8) 




53390.0 


52473- 


-54240 


129 


340 




1.48477368242(2) 


-0.017003(1) 




54400.0 


54262- 


-54821 


39 


306 


B1907-I-10 


3.52558388053(5) 


-0.03286(3) 




52578.0 


52470- 


-52686 


18 


161 




3.525581831083(6) 


-0.0328630(5) 




53300.0 


52705- 


-53689 


71 


391 




3.52557928013(1) 


-0.0328110(3) 


-2.02(4) 


54200.0 


53701- 


-54821 


88 


350 


B1913-I-10 


2.471906137401(5) 


-0.0931843(1) 


-0.071(9) 


53308.0 


52473- 


-54159 


118 


311 




2.47189777002(2) 


-0.093219(1) 




54348.0 


54180- 


-54831 


61 


306 


B2224-I-65 


1.465111203809(5) 


-0.0207375(4) 


-0.268(9) 


53880.0 


52487- 


-54265 


220 


922 




1.46511023680(2) 


-0.020737(2) 




54420.0 


54275- 


-54831 


71 


484 
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Figure 1. The 2005 December glitch of PSR B0144+59: (a) tim- 
ing residuals relative to the pre-glitch model, (b) variations of the 
frequency residual Av relative to the pre-glitch solution, (c) the 
variations of v. The dashed vertical line indicates the epoch of 
the glitch. 

3.1 PSR B0144+59 

This pulsar has a period of 196 ms and a small period 
derivative implying a relatively large characteristic age of 
12.1 Myr. Normally it is a very stable pulsar and has no 
previously detected glitch. Despite this, we detected a very 
small glitch at MJD ~ 53682, with a fractional frequency 
step of just 0.055 x 10~^. There are only four measured 
glitches with fractional size smaller than this. Also, PSR 
B0144-I-59 is the second oldest pulsar for which a glitch has 
been observed, with the oldest being the millisecond pulsar 
PSR B1821-24 which has a characteristic age of 29.9 Myr. 
Fig.[T]shows that after the glitch there is continuous increase 
of the frequency relative to the extrapolated pre-glitch (i.e., 
positive A!>p or decrease in that lasts for at least 600 d. 
This is unusual and opposite to the normal post-glitch decay. 
It is possible that this tiny glitch in PSR B0144-(-59 is just a 
fluctuation in the timing noise, but it does appear distinct in 
the Ajy plot. Similar decrease s in |z>| following a glitch were 
observed for PSR B19514 -32 (iFos tcr ct al.! ll990l ') and in the 
RRAT PSR J1819-1458 (|Lvne et al. 200c|). These may im- 
ply an increase in the moment of inertia or a decrease in 
the spin-down torque at the time of the glitch. For the PSR 
B0144-I-59 position fit given in Table [U whitening of the 
residuals required two frequency derivatives. 



3.2 PSR B0402+61 

This pulsar is also relatively old with Tc ~ 1.7 Myr. Fig. [2] 
shows a small glitch at MJD ~ 53041 (2004 February) with 
Ai/g/i/ - 0.6 X 10"^. In contrast to PSR B0144-I-59, there is 
an almost linear decrease of v relative to the pre-glitch value 
that lasts for at least 1000 days and overshoots the pre-glitch 
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Figure 2. The glitch of PSR B0402-I-61: (a) timing residuals rel- 
ative to the pre-glitch model, (b) variations of frequency residual 
At/ relative to the pre-glitch solution, (c) the variations of i>. 

extrapolatio n. This is similar to the behaviour seen in the 
Crab pulsar l|Lvne et al.lll99j ). The negative Ai>p is seen in 
Fig.[2](b) and is listed in Tabled 

3.3 PSR B0525+21 

This very long-period pulsar has had three small glitches re- 
ported. The first at MJD 42057 (January 1974) had a mag- 
nitude of /S.Vg/v — 1.2 X 10"' and an exponenti al decay with 
Q = 0.5 and timescale 150 days was measured (|Downslll98^ : 
emar fc Lvne|[l99^ ').[j anssen fc Stapper^ l|2006h reported 
a second glitch of similar magnitude at MJD 52284 and sug- 
gest that a very small glitch with Avg/u = 0.17 x 10~' 
may have occurred at MJD 53375. Our data, shown in 
Fig. [31 confirm the glitch at MJD 52284 (our central date 
is MJD 52280 but the dates are consistent within the un- 
certainties). The measured Q and Td values are similar to 
those for the 1974 glitch. Our data also show another small 
glitch of magnitude Aug/u = 0.43 x 10"' at MJD 53980. 
However, Fig. [3] shows no evidence for the small glitch 
(of ampli tude Au ~ 0.05 nHz) at MJD 53375 claimed by 
Q^anssen fc StappersI l|2006l ). Although this epoch coincides 
with a small gap in our data, a limit of 0.01 nHz can be 
placed on the ampl i tude of any glitch at this time. Since 
IJanssen fc Stapperi (|2006l ) did not show time or frequency 
residuals for this pulsar, it is not possible to comment on 
what led them to believe that there was a glitch at this 
time. We note that it is the smallest of the three glitches 
that they reported for this pulsar. 



3.4 PSR J0631+1036 

PSR J0631+1036, discovered bv IZepka et~aLl l| 19961 ). has a 
period of 288 ms and a very large period derivative of 105 x 
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Figure 3. Glitches of PSR B0525+21: (a) variations of frequency 
residual h.v relative to the pre-glitch solution, and (b) the fre- 
quency first derivative v. The long dashed lines are at the epochs 
of well-defined glitches; th e short dashed line is at th e epoch of 
the small glitch claimed by Ijanssen fc Stapperj ||2006| ) . 



10~ 



implying a small characteristic age ~ 44 kyr. As shown 



in Fig. m three normal glitches and one unusual event which 
may be a slow glitch were detected in our data. The first 
glitch at MJD 52852 has an amplitude of Aug/u ~ 19 x 
10~^ and shows a clear relaxation which is well described 
by an exponential decay model with Td ~ 120 d and Q ~ 
0.62 (Table [3]). A small second glitch was detected around 
MJD ~53229. Although it is not so obvious in Fig. 01 timing 
residuals around the time of the glitch given in Fig.[5]confirm 
its reality. Its relative magnitude is ~ 1.6 x 10"". 

Fig. 2] shows another event starting at MJD ~ 54129 
which may be related to the s low glitches dete c ted in PSR 
B1822-09 (|Shabanoval 1 19981 : IZou et all l2004l : IShabanoval 
l2007l ). At MJD ~ 54129 there is a sudden increase in z> (de- 
crease in \u\) leading to a gradual increase in i/ of ~ 23 nHz 
over about 130 d. During this time, v steadily decreased, 
over-shooting the pre-event value and returning v to the ex- 
trapolation of its behaviour before the event. Table [3] lists 
the maximum t^Vgjv of ~ 6.6 x 10~^. This behaviour is sim- 
ilar to that of the slow glitches observed in PSR B1822— 09 
IShabano va 2007) except that the rise in ii appears more 
abrupt and no overshoot is observed for PSR B1822— 09, 
leading to a more permanent increase in v relative to the 
pre-glitch solution. 

A fourth glitch occurred in June 2008, with an ampli- 
tude of Ai/g/i/ ~ 44 X 10~^. Although the data span after 
this event is short, it can be seen from the Fig. |3]that there 
is a jump in |!>| ~ 3 x 10~^* s~'^ at the time of the glitch fol- 
lowed by some recovery. It is too early to say if the recovery 
is exponential, similar to that for the glitch at MJD 52852. 



3.5 PSR J1705-3423 

This p ulsar, discovered in the Parkes Southern Pulsar Sur- 
vey bv lManchester et al.l l| 19961 ). has a period of 0.255 s and 
period derivative of 1.07 x 10~^^. A small glitch was detected 
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Figure 4. Glitches of PSR J0631-I-1036: (a) timing residuals rel- 
ative to the pre-glitch model, (b) variations of frequency residual 
At/ relative to the pre-glitch solution, and (c) the frequency first 
derivative v. 
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Figure 5. Timing residuals around the time of the small glitch 
in PSR J0631-I-1036. 



in Oct 2007. A plot of the time variation of Ai/ is given in 
Fig. E^. The dominant effect is a small change in v, with a 
magnitude of ~ 2.0 x 10~^ Hz. Fig. |61d shows no obvious 
change in v. 



3.6 PSR B1737 30 

PSR B1737— 30 is a young radio pulsar which has ex- 
hibited 20 glitches in 20 years, maki ng it one of the 
most frequently glitching pulsars known ('McKenna fc Lyng 
1990 ; Shc ma r fc L yncJ996. ;.Krawczvk et al,..2003; Zou et all 
2008i ). The relative glitch amplitudes range widely from a 
few parts in 10~^ to nearly 2 x 10~^. Significant increases 
in \i>\ at the time of the glitch are observed, at least for the 
larger glitches, and these tend to decay linearly with time. 
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Figure 6. Timing beliaviour of PSR B1705-3423: (a) variations 
of rotational frequency h.v relative to the pre-glitch value, and 
(b) variations of the frequency first derivative v. 



FIr. [7 ] shows the last (giant) ghtch observed by IZou et aP 
lj2008t ) and two further smaller glitches at MJD ~ 54447 
and ~ 54694. Following the large glitch there is a short- 
term (~ 100 d) exponential recovery in followed by an ap- 
proximately linear increase in v until the next glitch. This 
behaviour is very similar to that observed in the Vela pulsar 
jLvne et al.lll99tj ) as will be discussed further in [21 Follow- 
ing this large glitch, no further glitches were observed for 
1400 days, an unusually long "quiet" period for this pulsar. 
A smaU glitch with iS.Vg/v ~ 41 x 10"^ was detected at MJD 
~ 54447. Then, about 280 d later, we observed another small 
glitch with tSvgjv ~ 2.4 x 10~®. This brings the mean ghtch 
rate back close to one per year. 



3.7 PSR J1751-3323 

T his pulsar was discov ered in a Parkes multi-beam survey 
bv lKramer et al.] (|2003l ). In addition to its noisy timing be- 
haviour, two small glitches were detected around MJD 53004 
and 54435 with t^Vgjv ~ 2 x 10"^ and l^Vgjv ~ 3 x 10"'' 
respectively. Fig. [8] shows that no significant changes in i> 
can be associated with these glitches. 




500 1000 1500 

Days after MJD 52900 



2000 



Figure 7. Glitches of PSR B1737— 30: (a) variations of rotational 
frequency Ai^ relative to the pre-glitch value, (b) an expanded 
plot of Ai/ where the mean post-glitch value has been subtracted 
from the post-glitch data, and (c) variations of the frequency first 
derivative ii. 
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Figure 8. The glitches of PSR J1751-3323: (a) timing residuals, 

(b) frequency residuals t^v relative to the pre-glitch solution, and 

(c) the variations of v. 



Table 3. The glitch parameters. 
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Extrapolated Fitted h-^ 



Pulsar 
Name 


rLpocn 
IVIJJJ 


A.Ug/u 

(10 ) 


/I n— 3\ 
(10 j 


AUg/u 

(-t A— 9^ 
(lU ) 


AUg/v 

/I n— 3\ 
(10 ) 


AUp 
(10 s j 




w 


res 

(/iSj 


Data Range 


B0144+59 


53682(15)* 


0.056(3) 


-0.21(5) 


0.055(8) 


-0.22(8) 


0.0015(8) 


- 


- 


128 


52486- 


-54831 


B0402+61 


53041(6)* 


0.62(4) 


0.5(1) 


0.60(4) 


0.47(16) 


-0.007(2) 


- 


- 


407 


52469- 


-54517 


B0525+21 


52280(4)* 


1.53(5) 


0.9(1) 


1.6(2) 


1.1(1) 




0.44(5) 


650(50) 


620 


51547- 


-53975 




53980(12)* 


0.5(1) 


0.13(12) 


0.43(6) 


0.24(18) 








769 


53387- 


54517 


J0631+1036 


52852.0(2)* 


17.6(9) 


3.3(2) 


19.1(6) 


3.1(6) 




0.62(5) 


120(20) 


518 


52606- 


-53143 




53229(10) 


1.5(4) 


-0.08(5) 


1.6(7) 


0.3(4) 




- 


- 


249 


53051- 


-53387 




54129-54260 


6.6(4) 


-2.5(2) 




















54632.41(14)* 


44(2) 


5.7(5) 


44(1) 


4(2) 




0.13(2) 


40(15) 


257 


54437- 


-54750 


J1705-3423 


54384(10)* 


0.51(4) 


-0.1(5) 


0.52(5) 


0.0(6) 








660 


52486- 


-54832 


B1737— 30 


54447.41(4)* 


42(3) 


0.2(4) 


41.0(2) 


0.7(1) 








1811 


53101- 


-54696 




54694(3)* 


2.2(4) 


0.61(5) 


2.4(3) 


0.86(6) 








299 


54463- 


-54807 


J1751-3323 


53004(6)* 


2.3(3) 


-1.4(6) 


2.0(2) 


-2.1(4) 








883 


52496- 


-53433 




54435(8)* 


3(1) 


-1(3) 


3.0(4) 


1(3) 




— 


— 


3683 


53051- 


-54620 


B1758-23 


53309(18) 


493.2(4) 


0.26(5) 


494(1) 


0.19(3) 




0.009(2) 


1000(100) 


2896 


52494- 


-54606 


B1800— 21 


53444(16) 


3910(12) 


8.69(8) 


3914(2) 


10.0(4) 


—48(2) 


0.009(1) 


120(20) 


618 


53218- 


-53778 


B1809-173 


53105(2)* 


14.4(5) 


3(1) 


14.8(6) 


3.6(5) 




0.27(2) 


800(100) 


1080 


52529- 


-54394 




54367(13)* 


1.5(3) 


0.4(9) 


1.6(4) 


0.4(9) 








1070 


53450- 


-54764 


B1815-14 


52057(7) 


0.54(5) 


-0.7(4) 


0.54(4) 


-0.8(4) 








511 


51510- 


-52498 


B1822-09 


53719-53890 


7.2(3) 


-15.8(5) 




















54115.78(4)* 


122(2) 


-1.8(5) 


121.3(2) 


-1.9(3) 




— 


— 


935 


53742- 


-54395 


B1823-13 


53238.2(7)* 


3.4(3) 


0.07(1) 


3.5(3) 


0.068(6) 








584 


52518- 


-53730 




53734(5) 


2416(4) 


13.0(4) 


2417(2) 


14(5) 


-44(9) 


0.015(4) 


75(25) 


2668 


53509- 


-54080 


B1838-04 


53408(21) 


578.7(4) 


1.7(3) 


578.8(1) 


1.4(6) 


-0.21(1) 


0.00014(20) 


80(20) 


340 


52536- 


-54306 


J1853+0545 


53450(2)* 


1.49(5) 


1.6(5) 


1.46(8) 


3.5(7) 




0.22(5) 


250(30) 


333 


52494- 


-54131 


B1900+06 


54248(11)* 


0.33(3) 


0.04(9) 


0.33(3) 


0.04(9) 








334 


52473- 


-54821 


B1907+10 


52700(16)* 


0.27(7) 


0.5(7) 


0.27(13) 


0.5(10) 








201 


52470- 


-53198 




53700-54400 


1.52(5) 


-5.0(7) 


















B1913+10 


54162(1)* 


2.55(3) 


0.08(5) 


2.52(9) 


0.3(1) 








294 


52473- 


-54516 


B2224+65 


54266(14)* 


0.36(8) 


-0.8(4) 


0.39(7) 


-0.6(4) 








417 


53600- 


-54706 



* Glitch epoch determined by phase fit 
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Figure 9. The glitch of PSR B1758-23; (a) variations of fre- 
quency residual Av relative to the pre-glitch solution, (b) an 
expanded plot of Ai/ where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the frequency first 
derivative !>. 



3.8 PSR B1758 23 

This pulsar has a period of 415 ms and large period deriva- 
tive giving it a small characteristic age of 58 kyr. PSR 
B1758— 23 is distant and has a large dispersion measure (DM 
= 1074 cm~^ pc), resulting in a scattered pulse profile at 1.4 
GHz. Furthermore, the pulsar suffers frequent glitches, so it 
is difficult to obtain a position measurement from timing; 
the value quoted in Table [T] is from the VLA observations of 
iFrafi et al.l (|l993l ). 

Six glitches in this pulsar have been reported 
jKaspi et al.ll l993': She mar fc Lto3|1996I : IWang et al.ll2000l : 
iKrawczvk et al. ,2003. ), with A.Vg/u in the range (0.02 - 
0.34) xlO~®. Little or no change in v is observed, showing 



that glitches in this pulsar have no significant post-glitch re- 
covery. Our observations (Fig. [9| reveal a further glitch with 
/^Vg/u ~ 0.494 X 10"*^ at MJD ~ 53309, the largest seen for 
this pulsar. As for previous glitches in this pulsar, there is 
little or no post glitch decay (Fig.[S}3 and c), with the fitted 
value of Q ~ 0.009. 



3.9 PSR B1800 21 

Two giant glitches were detected f or PSR B1800— 21 in 
1990 and 1997 w ith /\ug/v ~ 10"® (|Shemar fc Lvndll996l : 
IWang et"aLll2000l '). while a small glitch with tWg jv ~ 5 x 
10"^ in 1996 was found by IKrawczvk et al.l (|2003[ ). Another 
giant glitch was detected in 2005 March (MJD ~ 53444) in 
our observations. As shown in Fig. 1101 the fractional jump 
in frequency is I^Vglv ~ 3.9 X 10"^. There is an exponen- 
tial recovery of a small part of the change in v and ii with 
a time constant of about 120 days, but the main observed 
recovery is linear. A similar linear increase in v following the 
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Figure 10. The glitch of PSR B1800-21: (a) variations of fre- 
quency residual Ai/ relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the variations of v. 

previous jump can be seen in Fig. [10] (c) . Once a gain, this is 
similar to the Vela pulsar post glitch behaviour (jLvne et al.l 
1996). Table [2] shows a large v both before and after the 
glitch corresponding to this linear recovery. 

3.10 PSR B1809-173 

PSR B1809— 173 suff'ered a modest jump of fractional size 
Ai/g/i/ ~ 14.3 X 10~^ and /Avjv ~ 3.3 x 10"^ in 2004 April 
(MJD ~ 53105). Observed frequency variations are shown 
in Fig. 111! An exponential decay model with Td ^ 800 d 
and Q ~ 0.27 is an adequate representation of the post- 
glitch relaxation. A second small glitch at Sep 2007 (MJD 
~ 54367) with /Avgjv ~ 1.6 x 10^** is also detected. 

3.11 PSR B1815 14 

Fig. [12] shows that a small glitch was detected in PSR 
B1815-14 around MJD 52057. Fits of a second-order poly- 
nomial to the pre-glitch data and cubic polynomial to the 
post-glitch data are given in Table [2] Both these and the 
phase-coherent glitch fit are consistent with a glitch of mag- 
nitude lAvglv ~ 0.54 X 10"® and lAvjv ~ -0.8 x 10"^. The 
post-glitch data show a quasi-periodic oscillation in resid- 
uals with period ~ 1200 d as illustrated in Fig. 1131 Simi- 
lar quasi-periodic oscillations in timin g residuals have bee n 
observed in very long data spans by iHobbs etHI Hooi). 
An increase in v (decrease in \v\) was observed from MJD 
~ 54000 lasting for ~ 200 d. This event, which is similar to 
the slow glitches observed in PSR B1822— 09, results in an 
increase in v, with [\v ~ 5 nHz. It is possible that similar 
but smaller fiuctuations in v at earlier times result in the 
quasi-periodicity in v and the timing residuals. 
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Figure 11. The glitches of PSR B1809-173: (a) variations of 
frequency residual Ai/ relative to the pre-glitch solution, (b) an 
expanded plot of At' where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the frequency first 
derivative !>. 
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Figure 12. The glitch of PSR B1815— 14: (a) timing residuals rel- 
ative to the pre-glitch solution, (b) variations of frequency resid- 
ual Ai^ relative to the pre-glitch solution, and (c) the variations 
of V. 
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Figure 13. Timing residuals for PSR B1815— 14 after the glitch 
with respect to the model given in Table [2] 



3.12 PSR B1822 09 

A total of seven glitches in PS R B1822— 09 have been re- 
ported in the previous 12 years (|Shabanovalll998l : IZou et al] 
12004 IShabanovallioOTi ). Five of these glitches are "slow", 
characterised by a relatively sharp rise in v followed by a 
decline to roughly its pre-glitch state over 100 - 200 d. This 
results in a gradual spin-up of the pulsar relative to its pre- 
glitch state to a new level at which it remains till the next 
event. A similar event was observed in PSR J0631-I-1036 
as reported abov e in § 13.41 Here we extend the analysis of 
IZou et al] (|2004l ) by another four years. Fig. [14] shows that 
a slow glitch started at MJD ~ 53720 and a "normal" glitch 
occurred a t MJD ^54116 . The former was also seen by 
IShabanoval (|20o3) around MJD 53740 near the end of her 
data set. The frequency first derivative z> shows a sudden 
jump of ~ 1.6 X 10~^^ at the start of the event. The spin 
frequency increases for about 200 d with respect to the pre- 
glitch fit until the increase in v decays away. 

For the larger and latest glitch, l^Vgjv ~ 0.121 x 10"''. 
Probably because of timing noise and the short post-glitch 
data span, no cl e ar po st-glitch exponential recovery can be 
seen. IZou et ahl l|2004h suggested that a small-magnitude 
slow glitch probably occurred near the end of their data 
span around MJD 52900. Although there is a small fiuctu- 
ation in v and v visible in Fig. [14] at this time, we consider 
that it is not clear enough to be identified as a slow glitch. 



3.13 PSR B1823 13 

PSR B1823— 13 is a young pulsar (~ 21.4 kyr) with a pe- 
riod of 10 1.4 ms. Two large glitches have been previously re- 
ported bv lShemar fc Lvnel ljl996l ) with l^Vgjv ~ 2.7 x lO"*^ 
and I^Vglv ~ 3.0 x 10"^ in 1986 and 1992 respectively. Both 
had a partial exponential recovery with timescale ~ 100 d 
followed by a linear increase in v. As shown in Fig. 1151 two 
new glitches were detected in the Nanshan timing. A small 
glitch with Ai/g/i/ ~ 3.5 x 10"^ occurred at MJD ~ 53238, 
and a large glitch at MJD ~ 53734 with /^Vgjv ~ 2.4 x 10"*^. 
There is no obvious relaxation after the small glitch, but af- 
ter the large glitch the relaxation is very similar to that of 
the earlier large glitches in this pulsar. Fitting for a single ex- 
ponential decay model with ^ 75 d and Q ~ 0.015 gives 
a relatively good fit to the post-glitch data between MJD 
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Figure 14. The glitches of PSR B1822-09: (a) variations of fre- 
quency residual Au relative to the pre-glitch solution, (b) an ex- 
panded plot of Au where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the variations of !>. 



Days after MJD 51 600 

Figure 15. The glitches of PSR B1823-13: (a) variations of fre- 
quency residual Au relative to the pre-glitch solution, (b) an ex- 
panded plot of Au where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the frequency first 
derivative u. 



53509 and 54080. Because of the small Q, there is effectively 
a permanent change in !>, Ai^p ~ -44 X lO"^'^ s"^ Fig. [T5l 
shows a linear increase in i> both before (through the small 
glitch) and after the larger glitch. Fitting to the linear part 
of the pre- and post-glitch data, despite the short data span 
after the glitch, we have v approximately 100 x 10"^^ 
and 65 x 10~^^ s~^ respectively. We will further discuss this 
in 31 

3.14 PSR B1838-04 

A large glitch with Aug/u ~ 0.58 x 10~^ was observed in 
PSR B1838-04 between MJD 53387 and 53429. This is the 
first glitch o bserved in this middle -aged pulsar which was 
discovered bv lClifton fc Lvnd (|l986t ). Fig. ll6l shows a hint of 
a short exponential decay, but the main effect of the glitch 
is what seems to be a permanent change in />. Fitting for 
this with a 80 d exponential decay gave a value of Aup of 
-0.21 X 10"^^ s-2 (TableEl). 

3.15 PSR J1853+0545 

We have detected the first kn own glitch for this recently dis- 
covered (|Kramer et al.ll2003f ) and relatively old pulsar (3.3 
Myr). The glitch is small with a magnitude of Avgjv ~ 
1.46x 10~^. Fig. ll7l shows a partial exponential recovery; fit- 
ting for this gives ~ 250 d, and fractional decay Q ~ 0.22. 
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Figure 16. The glitch of PSR B1838-04: (a) variations of fre- 
quency residual Av relative to the pre-glitch solution, (b) an ex- 
panded plot of Av where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the variations of v. 



3.16 PSR B1900+06 

This pulsar has been monitored at Nanshan since 2002 and 
a small glitch was detected in May 2007 (MJD ~ 54248). 
The fractional size is Avgjv - 0.33 x 10"^. Fig. [TSl shows 
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Figure 17. The glitch of PSR J1853+0545: (a) variations of fre- 
quency residual /^v relative to the pre-glitch solution, (b) an ex- 
panded plot of l\v where the mean post-glitch value has been 
subtracted from the post-glitch data, and (c) the frequency first 
derivative !>. 
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Figure 18. The glitch of PSR B1900-I-06; (a) variations of fre- 
quency residual /\v relative to the pre-glitch solution, (b) the 
frequency first derivative i>. 



that where may be some post-glitch recovery, but the data 
span is too short for further analysis. 



3.17 PSR B1907+10 

PSR B 1907-1- 10 has a period of 284 ms and is relatively old 
(tc ~ 1.7 Myr). Fig. [15] shows a tiny glitch at MJD ~ 52700 
with ^Vgjv ~ 0.27 X 10~^. Following that, there is evidence 
for a slow glitch between MJD ~ 53700 and ~ 54400. A 
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Figure 19. The glitches of PSR B1907+10: (a) timing residu- 
als relative to the pre-glitch solution, (b) frequency residual Ai/ 
relative to the pre-glitch solution, and (c) the variations of i/. 
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Figure 20. The glitch of PSR B1913+10: variations of (a) fre- 
quency residual Ai/ relative to the pre-glitch solution, (b) the 
frequency first derivative v. 



gradual but significant increase in frequency of about 5 nHz 
was observed with a corresponding perturbation in z>. Al- 
though this may be just a fluctuation in the timing noise, it 
does have the characteristics of the slow jumps observed in 
PSR B1822— 09 and some other pulsars. 



3.18 PSR B1913+10 

This pulsar has a period of 405 ms and characteristic age 
Tc ~ 0.4 Myr. Fig. l20l shows that a small glitch with Aug/i/ ~ 
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Figure 21. Glitches of PSR B2224+65: (a) timing residuals rel- 
ative to pre-glitch model, (b) variations of frequency residual Av 
relative to the pre-glitch solution, and (c) the variations of v. The 
sh ort dashed lines in th e Av plot mark the small glitches claimed 
bv ljanssen fc Stapperd ((2006) and the longer dashed line marks 
the glitch detected in our data. 



2.5 X 10"'' was observed at MJD ~ 54162. Although the 
post-glitch data span is small, there is clear evidence for an 
increase in |!>| indicating some post-glitch recovery. 



3.19 PSR B2224+65 

PSR B2224+65 is the hig h-velocity puls ar associated with 
the Guitar Nebula (Cord es et al.l Il993l ). It has a period 
of 0.682 s and a characteristic age of 1.1 Myr. A giant 
glitc h with At/p/t/ ^ 1 .7 x 10"'=' occurred in 1976 Octo- 
ber l|Backus et al.|[l982l ) and three tiny glitches (Afg/i/ < 
0.2 X 10"^ ) between 2000 Decem ber and 2005 March were re- 
ported bv ljanssen fc StappersI l(2006 l . The long dashed line 
in Fig. [21] marks another glitch detected in the Nanshan 
data. This event at MJD ~ 54266 has Avjv ~ 0.39 x 10"^. 
There is some evidence that |i>| decreased after the glitch, 
similar to PSR B0144-h59 (§|3ll. 

The short dashed lines mark the secon d and t hird small 
glitches claimed by Ijanssen fc Stapperi (|2006l ) at MJD 
52950 and 53434 for which the estimated Afg values were 
0.12 X 10"^ Hz and 0.28 x 10"*' Hz, respectively. (Nanshan 
observations commenced in 2002 August, after the time of 
the first small glitch.) Fig. ETJb) shows that the pulse fre- 
quency was increasing (relative to the pre-glitch solution) 
across the interval covered by these claimed glitches, but 
the increase appears smooth and no glitch larger than about 
0.1 X 10"^ Hz is evident in the data. As for PSR B0525-I-21, 
IJanssen fc StappersI (|2006l ) did not give phase or frequency 
residual plots for this pulsar, so it is not possible to assess 
the evidence for their claimed glitches. 




5 6 7 

log[Age (yr)] 

Figure 22. Distribution of known pulsars versus pulsar charac- 
teristic age with the distribution of glitching pulsars marked by 
the shaded area. 



4 DISCUSSION 

The study of pulsar glitches provides one of the few ways of 
probing the interior of neutron stars. Previous studies have 
shown that glitch properties, specifically characteristics of 
the post-glitch recovery, vary greatly from one pulsar to an- 
other and for different glitches in the same pulsar. Further- 
more, while inter-glitch intervals are relatively short (^1 yr) 
in some young pulsars, they are much longer in older pulsars, 
with many having no observed glitch over observing spans of 
30 years or more. It is therefore valuable to continue timing 
observations of a large sample of pulsars in order to bet- 
ter characterise glitch properties and hopefully to lead to a 
better understanding of the mechanisms involved and the 
underlying properties of neutron stars. 

This paper reports on glitches detected in continuing 
observations of a sample of 280 pulsars using the 25-m 
Nanshan telescope of Urumqi Observatory. We detected 29 
glitches in 19 pulsars, with this being the first detection of 
a glitch in 12 of these pulsars. With those already listed in 
the ATNF Pulsar Catalogue, this increases the number of 
pulsars with detected glitches to 63 and the total number 
of detected glitches to 199. Fig. l22l shows the distribution of 
glitching pulsars among all known pulsars as a function of 
characteristic age. About half of all known pulsars with ages 
less than 3 x 10'* yr have detected glitches, but for older pul- 
sars the fraction is much less, with no glitches being detected 
in pulsars older than 3 x 10^ yr. 

The range of fractional glitch amplitude ^Vgjv for 
glitches detected in our work is from 4.4x10"" 10 3.9x10"^^. 
Fig. 1231 shows the distribution of fractional glitch amplitude 
for all known glitche^. The detection threshold for a glitch 
depends upon the quality of the TO As, the frequency of the 
observations and the intrinsic timing noise. TOA uncertain- 
ties depend on the telescope and receiver parameters and the 
pulsar flux density and pulse width. For the Nanshan obser- 
vations, they are typically 100 - 300 /is. The pulsar timing 
project at Nanshan benefits from frequent and regular ob- 
servations on a large sample of pulsars. Although there are 



The three small glitches claimed bv ljanssen fc Stapper j ||2006| ') 
for which we see no evidence have been omitted from this and 
subsequent figures. 
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Figure 23. Distribution of relative glitch amplitudes of all known 
glitches. Those detected at Urumqi Observatory are marked by 
the shaded area. 
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Figure 24. Fractional glitch amplitude Aug/v versus character- 
istic ago. The dash lino marks the amplitude of Hug/u = 10 . 



some longer gaps due to instrumental or scheduling issues, 
observations of all pulsars are typically at 10-day intervals. 
The level of intrinsic quasi-continuous timing noise varies 
greatly from pulsar to pulsar. This and the overall sensitiv- 
ity of the observations to small glitches is best assessed from 
the Ai/ plots given for each pulsar. Typically, the sensitivity 
is in the range (0.1 — 1.0) x 10~^ Hz. Our work has revealed 
eight small glitches with relative size < 10~^. These stud- 
ies have extended the sample of small glitches and indicate 
that the frequency of occurrence of small glitches is com- 
parable to that o f "giant" glitches with Afg/i/ > 10~^ (cf., 
iLvne et al.|[200ol ). 

Fig. [24] shows glitch fractional amplitudes plotted 
against the pulsar age. Although there is large scatter, 
there is a clear tendency for the average /S.Ug/v to in- 
crease with age up to somewhere between 10* and 10^ years 
and then to decrease with increasing age. The increased 
detection of smal l glitches reported in this a nd other re- 
cent papers (Krawczvk et al.. :2003 ; Cognard fc Back er 2004|; 
Ijanssen fc StappersI 20061 ') has clarified this trend. Very 
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small glitches with l^Vg/v < 10 are only seen in pulsars 
wit h age greate r than 10^ yr. 

iLvne et al.l l|2000l) found that the jump in v at the time 



log[kl (s ^)] 

Figure 25. Plot of glitch activity parameter, Ag versus the abso- 
lute value of the spin-down rate v. The dash indicates equal values 
of these two quantities. Upper limits of Ag are given by assum- 
ing one glitch Hvg/v = 10~^ occurred for young pulsars with no 
detected glitch and < 10^ yr. Frequent glitching pulsars are 
marked specifically. 



of a glitch, A!>g, was related to the value of v, with the frac- 
tional change A!>g/i> being approximately constant. Table [3] 
shows that most of the measured values are about 10~^, but 
there is a scatter of up to two orders of magnitude in the 
ratio, even for different glitches in a given pulsar. 

The glitch activity parameter, Ag (Equation [1} is a 
way of quantifying the contribution of glitches to u. Since 
pulsars slow down, v is negative, but glitches are spin-ups 
and hence the contribution is positive. In some models, 
the fractional contribution is related to the fraction of the 
neutr on-star momen t of inertia which is in superfluid form 
(e.g., Ru derman et al. 1998 1. Fig. 1251 shows the dependence 
of activity parameter on |z>|. Again there is a large scat- 
ter, but there is a tendency for increasing activity with in- 
creasing |i>| with the ratio between them ranging between 
10~^ and 10~^. There seems to be a turn-over for the pul- 
sars with largest |i>j, for example the Crab pulsar and PSR 
J0537— 6910, having relatively smaller activity parameters. 
There is one pulsar, PSR J2301-I-5852 (the 7 s anomalous 
X-ray pulsar IE 2259-1-586 in the supernova remnant CTB 
109), where this ratio is close to 1.0. However, it should be 
not ed that only one glitch has been observed in this pul- 
sar l|Kaspi et al.ll200^ . so the derived activity parameter is 
effectively an upper limit. 

Glitches are generally attributed to the release of stress 
built up as a result of the steady spin-down of the pulsar. 
This str ess may be in the s olid crust, relieved by crustal star- 
quakes (|Bavm et al.lll96^ . or may be on pinned vortices in 
the superfluid i nterior , relieved by a collective unpinning of 
many vortices (A nderson fc It oh 1975). Both of these result 
in a sudden spin-up of the crust and hence a jump in the 
observed pulse frequency. In these models, one might expect 
a correlation of glitch size with the duration of the preceding 
or following inter-glitch interval, depending on the detailed 
mechanism. A strong correlation between the size of glitches 
and the length of the followin g inter-glitch interval wa s ob- 
served for PSRs J0537-691 jMiddleditch et al.ll2006l ) and 
B1642-03 (|Shabanovall2009l ). but such a clear correlation is 
not observed in other pulsars. 
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Figure 26. Fractional glitch amplitude versus length of the pre- 
ceding interglitch interval. 
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Figure 27. Fractional glitch amplitude versus length of the fol- 
lowing interglitch interval. 

Fig. [26] and [27] shows the dependence of the fractional 
glitch amplitude on the length of the preceding and fol- 
lowing interglitch intervals respectively for the most fre- 
quently glitching pulsars. Despite the naive expectation that 
glitches would be large after a long build-up period, there 
is very little correlation of glitch amplitude with length of 
the preceding interval. For example, the Crab pulsar and 
PSR B1737— 30 have a wide range of both parameters and 
there is essentially no correlation. Most glitches in the Vela 
pulsar are large and have long inter-glitch intervals. Two 
smaller glitches have been observed with shorter preced- 
ing intervals, but the preceding interval for the smallest 
glitch is more than four times the interval for the next 
smallest which is more than an order of magnitude larger. 
For PSR J0537-6910, there is no correlation and, for PSR 
B1338— 62, if anything there is a negative correlation. For 
po st-glitch intervals. Fig. [57] shows the correlations found 
bv lMiddleditch et al.l ( |2006 ) and .Shab anova (2009) for PSR 
J0537-6910 and PSR B1642-03 respectively, and a possi- 
ble steeper positive dependence for PSR B1338— 62. Again 
for the Crab and Vela pulsars and PSR B1737— 30 there is 
little correlation. 



Pulsars have a wide variety of post-glitch behaviours. 
Most show some form of relaxation toward the extrapolated 
pre-glitch solution after a glitch. However, there are some 
cases where the pulsar simply keeps spinning at the new 
rate with no change in i>. In terms of the exponential decay 
model (Equations [3] - [Zl , Ai/d ~ 0, Ai/g « Avp, Q ~ and 
Aug ~ 0. PSR B1758— 23 is a good example of a pulsar with 
large glitches having these characteristics. For small glitches, 
it is generally not possible to measure any decay since it is 
masked by random and systematic timing noise. 

In some pulsars there is an apparently permanent 
change in at the time of a glitch, for example, PSRs 
B0144-f59, B0402-I-61, B1838-04 and B2224+65. For PSR 
B0144+59 and possibly PSR B2224-f 65, the change is posi- 
tive which is unusual and not in accord with the exponential 
decay model. For the others, it is of course possible that the 
decay is exponential but with a time constant much longer 
than our post-glitch data span. 

Exponential recoveries are seen in many pulsars, but 
these generally have relatively short timescales, ~ 100 d, 
and low Q, that is, only a fraction of the initial frequency 
jump decays. Examples are the first observed glitch in PSR 
J0631-I-1036 and the large glitches in PSRs B1800-21 and 
B1823-13. For PSR B1838-04, there is a short exponential 
recovery preceding the long-term decay mentioned above. 

These small exponential recoveries are often asso- 
ciated with other post-glitch behaviours. For example, 
they precede linear increases in i> in PSRs B1737— 30, 
B1800-21 and B1823-13. The latter two pulsars also 
show an effectively permanent increase in jz/j at the time 
of the glitch, that is, the linear post-glitch recovery is 
approximately parallel to but below that existing before 
the observed glitch (see Figures [10] and llSp . Very simi- 
lar b ehaviour is observed in the Vela pulsar (|Lvne et al.l 
I1996D. and in PSRs B1046-58, B1706- 44 and B1727-33 
(jShemar fc Lvnelll996l : [Wang et al.ll2000l '). Table [2 gives the 
observed !> slopes, i.e., i), for the linear portion of the recov- 
ery. The slope in the Vela pulsar is an order of magnitude 
greater than for the other pulsars. For PSRs B1800— 21 and 
B1823— 13, multiple glitches with linear behaviour are ob- 
served, but the slopes are distinctly different for the different 
glitches, varying by about a factor of two. Similar slop e vari- 
ations are evident in the Vela data ijLvne et al.lll99^ '). The 
apparent braking indices 



corresponding to the linear increases in i> are also given in 
Table [4] Values range between 12 and 40, emphasising that 
these changes are due to internal dynamics in the neutron 
star, not to magnetic-dipole braking (for which n = 3). Al- 
though the Vela pulsar has the largest values, there is no 
clear correlation of either the slope or the braking index 
with pulse period or characteristic age. Table[3]also gives the 
maximum observed fractional glitch size, showing that all of 
these pulsars have giant glitches. However, it is worth noting 
that linear recoveries are not observed in PSRs B1338— 62, 
J1617-5055, B1757-24 and J202H-3651, despite these hav- 
ing similar characteristic ages and similar maximum glitch 
sizes. 

Fig. l28l shows the fraction of Aug which decays, Q, ver- 
sus characteristic age. Except for one event with Q ~ 0.1, the 
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Table 4. Rate of change of u for pulsars which show a hnear post-ghtch recovery. 



Pulsar P Age Max i^Vgjv v n Note 



Name 
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(kyr) 


(10-6) 


(10-24 s-3) 












B0833- 
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0.089 


11.3 


3.06(7) 


873(10) 


40.3(5) 


MJD 


50024- 


-50364 


("Wans et al. 2000) 


B1046- 


-58 


0.124 


20.3 


2.995(7) 


72(10) 


15(2) 


MJD 


47909- 


-48929 


(Wane et al. 2000) 


B1706- 
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124(1) 
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MJD 
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(Wane et al. 2000) 
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26 
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(Shemar & Lvne 1996) 
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Figure 28. The fraction parameter Q in which a glitch l\Vg 
decays, versus characteristic age. 
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Figure 29. Glitch decay time-scale versus pulsar characteristic 
age. 



Crab pulsar usually has a large Q between 0.7 and 1.0. Al- 
though there is some tendency for decreasing Q with increas- 
ing age (Wang et al. 2000), several old pulsars, for example, 
PSRs B0525-f21, B1642-03, B1809-173 and J1853-h0545 
show some evidence for exponential relaxatio n. Fig . [29 



shows that, with the addition of PSR B1642— 03 (|Shabanovi 
l2009l ). there is a clearer correlation of decay ti mescale with 
age than was previously evident ((Wan g et al. 2000). This 
suggests a conne ction with the int ernal temperature of the 
neutron star (cf. lAlpar et al.|[l989l ). 

Two new examples of the "slow glitch" phenomenon, 
in PSRs J0631-fl036, and B1907+10, have been identi- 
fied in this paper. Previously, slow glitches have only been 
identified in one pulsar, PSR B1822 -09 (Shabanova 199i; 
Shabanova fc Uramall2000l: IZou et al.|[2004i : IShabanoval2005l . 



2007t ). Although it is possible to interpret these events as 
fluctuations in timing noise, they appear to be distinct and 
repeatable with similar characteristics. Hence it is reason- 
able to consider them as a separate category. PSR B1822— 09 
has had five slow glitches between 1995 and 2004 and 
one more in 2006 reported in this work. The new slow 
glitches detected in PSRs J0631-I-1036 and B1907-f 10 sig- 
nificantly increase the sample of slow glitches. The event 
in PSR J0631-I-1036 is similar in size and characteristics to 
those in PSR B1822— 09, with an increment in v of about 
2 X lO"^'^ s-2 lasting ~ 100 d. However it differs in that 
the decline in v overshot the pre-event value, resulting in 
an effective spin-down of the period to approximately its 
extrapolated value. In PSR B1822— 09, the spin-ups persist 
until the next event. 

The results presented in this paper significantly increase 
the sample of known pulsar glitches. They further illus- 
trate the complex nature of glitches, with a great diversity 
of glitch sizes and modes of post-glitch recovery being ob- 
served, both in individual pulsars and across different pul- 
sars. They also strengthen the identification of slow glitches 
as a distinct phenomenon. While there are several different 
models for glitches, as yet none of these models satisfactorily 
account for this diversity and so the promise that glitches 
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have as probes of the interior of neutron stars has not yet 
been fully realised. Hopefully, these and future results will 
act as a stimulus for a better understanding of the mech- 
anisms involved and their implications for the physics of 
neutron-star interiors. 
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